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A STUDY OF THE MIXING OF HYDROGEN INJECTED 


NORMAL TO A SUPERSONIC AIRSTREAM* 

By R. Clayton Rogers 
Langley Research Center 

SUMMARY 

An investigation has been conducted on the penetration and mixing of hydrogen 
injected normal to a supersonic airstream with the ratio of jet dynamic pressure to 
free-stream dynamic pressure varied from 0.5 to 1.5. Hydrogen gas was injected at 
sonic velocity from a circular injector in a flat plate into a Mach 4.03 airstream at a 
stagnation temperature of 300 K and stagnation pressures of 13.6 and 20.4 atmospheres 
(1 atmosphere = 101.325 kN/m2). Corresponding Reynolds numbers per meter were 
6.19 x 10? and 9.28 x 10^, respectively, and resulted in a turbulent boundary-layer thick- 
ness of 2.70 injector diameters at the injector station. Measurements of hydrogen vol- 
ume fraction and pitot and static pressures were made at downstream stations of 7, 30, 
60, 120, and 200 injector diameters and yielded mass fraction, Mach number, stagnation 
pressure, velocity, and mass flux profiles. Data correlations from the literature under- 
predicted the measured jet penetration; this underprediction is thought to be due, in part, 
to the relatively thick boundary layer of this investigation. At all downstream stations 
the penetration was found to be proportional to the 0.3 power of the dynamic-pressure 
ratio and the decay of the maximum mass concentration (to values of 0.1) was inversely 
proportional to the square root of the downstream distance. At any particular station 
the maximum mass concentration was proportional to the 0.54 power of the dynamic- 
pressure ratio. Nondimensional concentration profiles for the vertical surveys were 
approximated by Gaussian-type functions and showed similarity at stations equal to or 
greater than 60 injector diameters. 


INTRODUCTION 

Advanced hypersonic vehicles, such as a reusable launch vehicle with an airbreathing 
propulsion system for the first stage, are currently under investigation. (See ref. 1.) The 


*The information presented herein was included in a thesis entitled "The Penetration 
and Mixing of a Sonic Hydrogen Jet Injected Normal to a Mach 4 Airstream" offered in 
partial fulfillment of the requirement for the degree of Master of Science in Aerospace 
Engineering, Virginia Polytechnic Institute, Blacksburg, Virginia, March 1970. 


proposed propulsion system is a hydrogen-fueled supersonic combustion ramjet with oper- 
ation to a free -stream Mach number of 12, and a corresponding combustor entrance Mach 
number of approximately 4. The accurate analysis and design of the combustor requires 
a knowledge of the fuel-air mixing characteristics and a means of injecting the fuel so 
that a nearly uniform fuel distribution and a short mixing length are obtained without pro- 
ducing significant thrust penalties. These requirements suggest that parallel injection 
from multiple points in the stream (refs. 2 and 3) and normal injection from several 
points on the walls (ref. 4) will both be included in the optimized design. 

Normal sonic injection of various gases and gaseous mixtures issuing from a dis- 
crete circular hole in the flat plate into a supersonic airstream has been investigated and 
reported in references 4 to 10. These data are generally for conditions corresponding to 
values of the ratio of jet dynamic pressure to free- stream dynamic pressure greater than 
1.3 and indicate that for sonic injection, shorter mixing lengths occur at lower values of 
the ratio of jet dynamic pressure to free- stream dynamic pressure. Analytical methods 
presented in references 8 to 10 for predicting the penetration of the jet and the jet trajec- 
tory in the unconfined supersonic mainstream have been developed from empirical and 
semiempirical data correlations. Generally, these correlations do not consider any 
effect of the mainstream boundary-layer thickness. 

The present investigation was conducted to provide detailed information on the pen- 
etration and mixing of hydrogen injected normal to a supersonic airstream with a thick 
boundary layer for values of the ratio of jet dynamic pressure to free-stream dynamic 
pressure from 0.5 to 1.5. These data are needed to aid in constructing analytical methods 
that yield reliable predictions of the mixing process in supersonic combustors. 

The hydrogen gas was injected from a 0.102-centimeter-diameter sonic nozzle 
perpendicular to the surface of a flat plate mounted in a 23-centimeter-square tunnel 
test section. The tests were conducted at a free-stream Mach number of 4.03, stag- 
nation temperature of 300 K, and stagnation pressures of 13.6 and 20.4 atmospheres 
and resulted in Reynolds numbers per meter of 6.19 x 10^ and 9.28 X 10^, respectively 
(1 atmosphere = 101.325 kN/m 2 ). Boundary -layer thickness on the flat plate at the injec- 
tor station was 2.7 injector diameters. Ratios of jet dynamic pressure to free-stream 
dynamic pressure of 0.5, 0.75, 1.0, and 1.5 were investigated. Measurements of hydrogen 
volume fraction, pitot pressure, and static pressure were obtained by vertical and hori- 
zontal surveys of the flow field at downstream stations of 7, 30, 60, 120, and 200 injector 
diameters. The present data are compared with the previously obtained data correlations. 

SYMBOLS 


A 

bl 


2 


streamtube cross-sectional area, meters 2 
exponent defined in equation (7) 


exponent defined in equation (8) 
drag coefficient 

injector nozzle exit diameter, meters 
equivalent jet exit diameter, dK 1 /2 ) meters 
fuel-air mass ratio 

distance measured along center line of emerging jet, meters 
injector nozzle discharge coefficient 
Mach number 


mass flow rate, kilograms/ second 
Reynolds number 

absolute pressure, newtons/meter^ or atmospheres 

ratio of jet dynamic pressure to free-stream dynamic pressure, 

(p v 2 )j/(p v2 ). 

absolute temperature, kelvin 
velocity, meters/second 
longitudinal coordinate 
potential core length, meters 
lateral coordinate 
vertical coordinate 
hydrogen mass fraction 


air mass flux parameter, (pV) x (l - a), 


kilograms 
meter^- second 


boundary-layer thickness, meters 

boundary -layer displacement thickness, meters 


slope of emerging jet center line measured from horizontal, degrees 


6i initial slope of injector, degrees 

X ratio of jet mass flux to free-stream mass flux, (pV)jy/(pV) 

v hydrogen volume fraction 

£ hydrogen mass flow rate parameter, a(pV) x j(pV)j 

, p mass density, kilograms/meter 3 

Subscripts: 


°° ‘free-stream conditions 

0 conditions at edge of mixing region where v = 0.005 

1 conditions in undisturbed flow upstream of injector 

5 conditions at which the mass concentration is half maximum 

j jet conditions 

t stagnation conditions 

x survey point 

max maximum value 

ref reference value 

a conditions at which concentration is maximum 

A bar over a symbol denotes an average quantity. 


MODEL AND FACILITY 


Test Apparatus and Model 

' A sketch of the stainless-steel rectangular flat plate used in the experiments is 
shown in figure 1. A 0.1 02 -centimeter -diameter sonic nozzle was flush-mounted perpen- 
dicular to the plate surface 18.6 centimeters from the plate leading edge. Details of the 
nozzle are shown in figure 1. The plate leading edge was a 2° wedge tapering to a cylin- 
drical leading edge of approximately 0.013 -centimeter thickness, followed by a 10° wedge. 
The plate spanned the 23-centimeter by 23-centimeter test section of a continuous-flow 
supersonic tunnel. Additional information on the tunnel may be found in reference 4. 
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Boundary -layer profiles at the injector station with stagnation temperatures of 
300 K and stagnation pressures of 13.6 and 20.4 atmospheres (unit Reynolds numbers of 
6.19 x 10 7 and 9.28 x 10 7 per meter) are presented in figure 2. Also shown are theoreti- 
cal results, computed from the method of reference 11 for the 20.4 atmosphere condition, 
which agree well with the data. The boundary -layer thickness based on an edge velocity 
of 99 percent of the free -stream value was 2.70 injector exit diameters for both free- 
stream stagnation pressures. 


Injector Flow 

A schematic of the hydrogen gas supply and control system is presented in fig- 
ure 3. The system was equipped with pressure-controlled valves and purged with nitro- 
gen. Total temperature of the jet was measured by using a standard iron-constantan 
thermocouple inserted in a filter between the orifice meter and the injector. Jet total 
pressure was measured near the injector by means of a wall static-pressure orifice 
mounted in the 0.476-centimeter-diameter injector supply tube. Calculations indicated 
that the measured pressure would be within 99 percent of the jet total pressure. A tube 
attached to the injector supply line supplied 100-percent hydrogen samples for full-scale 
chromatograph readings. During a survey, gas samples of the hydrogen-air mixture 
were taken through the pitot probe with the aid of a vacuum pump at mass flow rates up 
to 4.17 mg/sec. The sample flow to the chromatograph was metered to 0.0556 mg/sec 
and the remaining flow bypassed and discharged to the atmosphere. The apparatus was 
operated over a jet total-pressure range of 2 to 4 atmospheres corresponding to ratios 
of jet dynamic pressure to free-stream dynamic pressure of 0.5 to 1.5. Test conditions 
are presented in the following table: 


Test 

condition 


Pt,oo> 

atm 

P M’ 

atm 

g/sec 

N Re,°°> 

m-1 

6/d 

1 

0.50 

20.4 

1.986 

0.0820 

9.28 X 10 7 

2.70 

2 

.75 

20.4 

2.980 

.1230 

9.28 

2.70 

3 

1.00 

13.6 

2.647 

.1094 

6.19 

2.70 

4 

1.50 

13.6 

3.960 

i 

.1641 

6.19 

2.70 


Instrumentation 

» 

Gas analyzer .- The volumetric concentration of hydrogen in the gas samples was 
measured by a process gas chromatograph. (See refs. 12 and 13.) At the beginning of 
a 1 -minute cycle, part of the sample gas was isolated and forced by the carrier gas, 
nitrogen, through a molecular sieve and a column consisting of a length of stainless -steel 
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tubing packed with silica gel. This procedure provided a qualitative identification of each 
component since each will process through the column at a predictable rate. The quan- 
tity of each component was determined by four thermal conductivity detectors of which 
two are always exposed to the carrier gas. The unbalance of the detector bridge pro- 
vided a voltage output proportional to the cooling effect and hence was a measure of the 
concentration of the separated sample components relative to the carrier gas. The volt- 
age output was recorded by a pen deflection on a strip chart. Readout controls were 
adjusted so. that only the hydrogen concentration was detected. 

Probe description .- The gas -sampling pitot probe and the static probe are shown in 
figure 4. The pitot -sampling probe is a boundary- layer survey type with the probe tip 
mounted in a 7.94-millimeter-diameter supporting tube offset to allow for actuator rod 
clearance. The actuator mechanism provided for probe movement for vertical traversing 
and yaw in the horizontal plane. The static-pressure probe was of similar design with a 
cone angle of 28° and four 0.203-millimeter orifices located at 14 probe diameters from 
the tip. 

Flow measurement .- The mass flow of the injected gas was measured with a 
0.318-centimeter-diameter, sharp-edged, corner-tap orifice meter. The static temper- 
ature at the meter was assumed to be the same as the jet stagnation temperature. Hydro- 
gen mass flow rate through the orifice meter was calculated from an equation derived from 
a hydrogen- corrected air calibration of the orifice meter. The discharge coefficient of the 
injector nozzle, based on orifice meter measurements, normally ranged from 0.73 to 0.78 
with an average value of 0.76. The sample flow rate to the chromatograph and the bypass 
flow rate were measured by thermoconductivity mass flow rate meters. All pressures 
except tunnel-wall static pressures were measured with strain -gage -type transducers 
and recorded on automatic balance potentiometers. The tunnel-wall static pressures 
were read on mercury manometers and recorded periodically during each test run. 

Survey Procedures 

At each of the downstream stations, one vertical and three horizontal surveys were 
made of the flow field. The vertical survey was made along the jet center line stepwise 
from the plate surface outward until a zero hydrogen concentration was obtained. Hori- 
zontal surveys were then made at points above the plate corresponding to maximum and 
half-maximum concentration and at a point midway between the plate surface and the point 
of maximum concentration. For each horizontal survey, the edge of the mixing region was 
located, and a stepwise survey was made across the flow field from this point until a hydro- 
gen volume fraction of zero was obtained. At each point in the survey, a gas sample and a 
pitot-pressure measurement were taken. 
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Data Reduction and Accuracy 

The pressure and concentration data at each survey point were reduced to yield val- 
ues of mass fraction, mixture molecular weight, Mach number, mixture total and static 
temperature, velocity, mixture and air mass flux, and the hydrogen mass flow parameter. 
The molecular weight of the hydrogen-air mixture was computed by assuming molecular 
weights of the components of 2.016 and 29.0, respectively. Mixture total temperature was 
obtained from the mixture mass averaged total enthalpy computed from the measured total 
temperatures of the jet and free stream. The mixture was assumed to be a perfect gas 
and the values of Mach number, static temperature, and velocity computed using the equa- 
tions for one -dimensional isentropic flow presented in reference 14. Local density of the 
mixture, used to compute the mass flux parameters, was calculated by using the ideal gas 
quation of state with the universal gas constant equal to 8.31 J/mol-K (1.986 cal/mol-K). 

During the course of the tests, it was determined that the center line of the hydrogen- 
air flow field was not always coincident with the tunnel center line. This difference is 
believed to be a result of the small scale of the injector and probe tip and slight asymme- 
tries of the tunnel flow. The maximum concentration was, therefore, sometimes obtained 
at a point to either side of the vertical survey location. The lateral location of the maxi- 
mum concentration was taken as the center line of the mixing flow and the vertical survey 
was considered to have been made at a point slightly off center. In most cases, the dis- 
tance between the tunnel and flow-field center lines (y/d) a was less than one injector 
diameter. Probe position accuracy of the actuator mechanism for the vertical surveys 
was ±0.127 mm, the same as the probe tip height. This value corresponds to a possible 
error in the vertical probe position of ±0.125 injector diameter. Accuracy in the yaw 
mode was ±0.10° over a yaw angle range of less than ±10°. At the widest survey location, 
the change in the x-position was less than two injector diameters. The gas chromatograph 
was periodically calibrated with 100-percent hydrogen from the supply line and the repeat- 
ability checked to a variation of less than one-half-percent full scale. The calibration of 
the instrument was nearly linear so that the error in any sample would be ±0.005 volume 
fraction. 


RESULTS AND DISCUSSION 

* 

Flow -Field Structure 

The general structure of the flow field resulting from the normal injection is pre-* 
sented in figure 5. The data shown in figure 5 are profiles of hydrogen mass fraction 
taken on the tunnel center line and show trajectories of various concentrations. The 
hydrogen jet is quickly turned downstream by the free stream and mixes rapidly near 
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the jet, the maximum concentration decreasing to about 12-percent mass fraction in 
7 injector diameters. Farther downstream, the mixing is slower, the maximum concen- 
tration decreasing to a mass fraction of 4 to 5 percent at an x/d of 60. The bow shock, 
determined from schlieren photographs, was not appreciably affected by a change in q r 
and was essentially a Mach line downstream of an x/d of 20. Details in the vicinity of 
the injector were not clear from the present tests; however, flow details in this region 
are presented in references 6,7, and 9 for thinner boundary layers and larger q r than 
the present, investigation. 

The extent of the separation depends on the boundary- layer thickness relative to the 
injector diameter and the degree of underexpansion of the jet. From physical considera- 
tions of the injection disturbance, it may be reasoned that with a constant jet diameter, 
injection into a thick boundary layer will produce a greater absolute penetration near the 
jet, but the jet will be turned downstream somewhat before encountering the high velocity 
free stream and a weaker bow shock in the free stream will result. Therefore, the injec- 
tor effective back pressure would be less. 

For the free-stream conditions of figure 5, the minimum jet condition for sonic 
injection was a q r of approximately 0.45. This value is equivalent to a jet exit static 
pressure and thus an effective back pressure of 0.63 atmosphere or approximately 40 per- 
cent of the free-stream pitot pressure. For these tests most of the jet penetration into 
the airstream occurs within 7 injector diameters, and the mixing region remains almost 
entirely embedded in the boundary layer. Variation of the static pressure vertically 
across the mixing region was less than 2 percent. 

Jet Penetration 

Although jet penetration has been discussed considerably in the literature, there are 
various definitions of penetration. As used herein, the term "penetration" is referred to 
the edge of the mixing region in the vertical center-line plane, where v is 0.005, and is 
denoted as (z/d) 0 ; the height at which the concentration is maximum is referred to as 
"penetration to a max " and is denoted as (z/d) a . Figures 6 and 7 compare the present 
penetration data with various data correlations from the literature and illustrate the effect 
of ( the dynamic -pressure ratio at an x/d of 7 and the maximum concentration and pene- 
tration trajectories for a q r of 1.0. 

t Correlations are presented in reference 9 for the penetration and a max trajec- 
tories and are given here for normal sonic injection of hydrogen in a Mach 4.03 free 
stream: 


(D 0 ■ 3 --(« r H5P 66 


( 1 ) 
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( 2 ) 



3.45 


( q r) 


0.533/x 
Vd 


-0.259 


These equations were derived for data at x/d equal to or greater than 7. At an x/d 
of 7 (fig. 6) , equation (1) underpredicts the penetration by about 20 percent and has a slope 
of 0.5 compared with 0.3 for a fairing of the experimental data. For a q r of 1.0 (fig. 7), 
the penetration trajectory predicted by equation (1) is as much as 45 percent lower than 
the data at an x/d of 200. Based on measurement accuracy, the probable error of the 
data points has been estimated to be less than 0.2 injector diameter or approximately 
3 percent. Equation (2), in figure 6, shows a reasonable agreement with the data and has 
a slope of 0.533 compared with 0.6 for a fairing of the data points. The a max trajec- 
tory given by equation (2), for a q r of 1.0 (fig. 7), approaches the plate surface at large 
values of x/d whereas the data diverge. 

Another correlation for the penetration of a normal jet, reported in reference 10 
and given here for sonic injection into a Mach 4.03 free stream, is 



0.392/x\0.134 


3.40(q r )°-^(§)' 


(3) 


Equation (3) was derived from data taken at free-stream Mach numbers of 1.6 and 3.0 
over an x/d range of 14 to 167. Equation (3) is presented in figures 6 and 7 and gives 
a better prediction of the penetration than equation (1). 

Also shown in figures 6 and 7 is the jet center-line trajectory ^ max trajectory) in 
the near field calculated from the method of reference 9 in which the jet is considered as 
being composed of cylindrical elements of length and the aerodynamic drag on each 

element is computed from empirical equations. The equation from reference 9 is 


where 



de 

C D (0) sin^fl 


(o.22 j* + 2.25) 4 


2.25 4 


(4) 


C D (0) = 1.2 + (m b sin e) 7/2 

(o £ Moo sin 0 = 

/ \ - 3 

/ 

C d ( 6>) = 1.06 + 1.14(Moo sin 6 ) 

(M^ sin 6 s 


Equation (4) underpredicts the effect of q r on the penetration at an x/d of 7. The 
effect of boundary layer is not accounted for in equation (4); however, calculations made 
by using the mass-averaged boundary- layer conditions rather than free-stream conditions 
made no significant change in the trajectory or the effect of q r> 
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Penetration Trajectories 


Trajectories of maximum concentration, half-maximum concentration, and the 
penetration height are correlated with the ratio of jet dynamic pressure to free- stream 
dynamic pressure and are presented in figures 8(a), 8(b), and 8(c), respectively. During 
examination of the data, it was found that the maximum concentration trajectories had a 

minimum value that occurred farther downstream at the higher values of q r . A factor 

1 fi * 

q r “ • applied to the x/d coordinate produced the family of curves presented in fig- 
ure 8(a). The turning of the maximum concentration trajectory beyond the parallel (that 
is, ^90°) with the plate surface, as shown by the initial decrease in (z/d) a , is thought 
to be more evident than in previous tests because of the thick boundary layer. Minimum 
values of (z/d) a occurred at values of x/d ranging from 9 to 45 for values of q r 
between 0.5 and 1.5, respectively. Examination of the trajectories of the point at which 
the mass concentration is one-half the maximum led to a correlating factor of q r “®'® 
which, when applied to the x/d coordinate, produced the family of curves similar to 
those for the maximum concentration, presented in figure 8(b). At all x/d stations, 
the penetration increased in proportion to the 0.3 power of the dynamic-pressure ratio. 
The correlated trajectories are presented in figure 8(c). For x/d less than 120, the 
data may be represented by the following equation: 




0.300/x\ 0 - 143 


(!) 


(5) 


This equation is similar to equation (3); the differences in the coefficient and exponents 
between equations (3) and (5) are probably due in part to the different injected gases and 
boundary -layer thicknesses. 


Decay of Maximum Concentration 

The concept of a potential core length has been used in reference 5 to nondimension- 
alize the longitudinal coordinate and to obtain a correlation of the maximum concentration 
decay with downstream distance. The potential core length Xq is defined as the distance 
along the center line from the injector to the downstream station at which the concentra- 
tion of the injected gas first decreases from 100 percent. For the present data, values of 
the potential core length, ranging between 0.14 and 0.25 injector diameter, were obtained 
by extrapolating the curves in figure 9(a) to a concentration of 100 percent and, therefore, 
must be considered as a correlating parameter rather than a true value of potential core 
length. Figure 9(b) presents the correlated maximum concentration decay data of the 
present investigation and includes the data for hydrogen from reference 5. The average 
deviation of the data from the correlated curve is less than 10 percent; the maximum 
deviation is 30 percent at an x/d of 7. The slope of the correlated curve for x/x 0 
less than 200 is approximately -0.5. Coaxial data from reference 3 indicated a maximum 
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concentration decay inversely proportional to x^. Reference 5 presented a correlation 
for the potential core length for normal injection as a function of the molecular weight and 
the ratio of specific heats of the injected gas, and the ratio of jet mass flux to free-stream 
mass flux. For hydrogen injection, the correlation reduces to 


^ = 0.317X 0,54 (6) 

d i 

and is compared with the present data in figure 10. Note that in equation (6), x Q is 
divided by the effective jet diameter dj which is the injector diameter corrected by the 
injector discharge coefficient. 


Profile Data 

Nondimensional profiles of concentration, velocity, and total pressure for the verti- 
cal survey and the horizontal survey through the point of maximum concentration are pre- 
sented in figures 11 to 13 and figures 14 to 16, respectively. The vertical concentration 
profiles in figure 11 are nondimensionalized and the origin of the coordinate system 
shifted to (z/d) re f. Note that the vertical coordinate near the wall is nondimensional- 
ized by a parameter different from that used above the point of peak concentration. The 
shape of the upper part of the profiles was not affected by q r and could be represented 
by a Gaussian- type function of the form 





where bj , the vertical profile shape index, is selected to give a reasonable fit with the 
data. At x/d stations, downstream of 30, the value of bj is constant at a value of 
2.70; thus, it is suggested that the flow field has become fully developed. It has already 
been noted that the maximum concentration trajectory in figure 8 has a minimum near a 
value of of 30. 

The section of the profiles at negative values of the vertical coordinate shows no 
systematic effect of q r ; however, the variation of concentration across the profiles is, 
in general, less at lower q r . 

Nondimensional velocity profiles, in comparison with the undisturbed boundary- 
layer velocity at the injector station, are presented in figure 12. For each data profile, 
the value V 0 is the velocity at the edge of the mixing region at a height (z/d) Q above 
the plate. Near the injector the peak velocity increases with increasing q r because of 
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the direct increase of the mass flow rate of injected gas with q r for constant free- 
stream conditions and a constant value of Vj. The initial ratio of jet velocity to free- 
stream velocity is approximately 2. At values of x/d downstream of 30, the shape of 
the velocity profiles approaches that of the boundary layer. 

Nondimensional total-pressure profiles in the vertical plane are presented in fig- 
ure 13 in comparison with the boundary-layer total-pressure profile. The total pressure 
at the edge of the mixing region generally increases with increasing x/d because of the 
addition to the mixing region of free- stream air that has passed through a weaker section 
of the bow shock. All the profiles exhibit approximately the same total-pressure loss 
due to the loss in momentum required to turn and accelerate the injected hydrogen. This 
region of low total pressure extends over approximately 60 percent of the vertical height 
of the flow field and is most severe at the upstream stations. 

Nondimensional horizontal concentration profiles are presented in figure 14 with 
the lateral coordinate nondimensionalized by the average distance from the center line to 
each of the lateral edges of the mixing region (y/d) Q . Values of (y/d) Q have a some- 
what random variation of roughly 5 to 8 from an x/d of 7 to 200 with a ±20-percent 
deviation at a given station. It was found that the profile data could be represented by a 
Gaussian-type function of the form 



where \>2 is the horizontal profile slope index and takes on a value required to obtain a 
reasonable fit with the data. 

Velocity profiles in the horizontal plane are presented in figure 15 and indicate a 
nearly uniform two-dimensional flow field downstream of about 60 injector diameters for 
q r of 1.0. Corresponding total-pressure profiles, presented in figure 16, show a large 
variation in total pressure across the mixing region because of the changing density of 
the mixture. The pressure level near the center line never recovers from the injector 
disturbance and is less than 8 percent of free- stream total pressure. The static- 
pressure variation across the mixing region is less than 2 percent. 

Flow -Field Contours 

An indication of the overall accuracy of the profile data was obtained by comparing 
the integrated hydrogen mass flow rate with the metered hydrogen mass flow supplied to 
the jet. In terms of £, the ratio of integrated mass flow to jet mass flow m int /mj is 
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( 9 ) 


”Hnt _ 4 ^max A o f 1 ' 0 | J A\ 

riij 7rdj^ ^0 ^max \Ao/ 

and is presented in figure 17. Agreement of the integrated and metered mass flows 
improves as x/ d increases and q r decreases. This improvement is probably due 
to the smaller gradients in the concentration and velocity - which are associated with 
the local turbulence level - at the downstream stations and lower values of q r . The 
characteristics of a binary gas flow field that affect the accuracy of concentration mea- 
surements are discussed in reference 3. For flow fields of this nature, differences 
between integrated and measured mass flow rates of 20 percent are considered to be 
typical. 

The contours of hydrogen mass fraction, presented in figure 18, represent cross 
sections of the flow field in the YZ -plane and are bounded by the zero concentration con- 
tour. At x/d stations between 30 and 200 and for all values of q r , the part of the con- 
tours above (z/d) a may be represented by semicircles centered at (z/d) a . This 
representation suggests that the flow field above (z/d) a is similar to coaxial mixing 
downstream of x/d of 30. 

The air mass flow rate contours in the nondimensional form £/|3 max and con- 
tained within the zero concentration contour are presented in figure 19. The mass flow 
rate of air contained within the mixing region was obtained by evaluating the integral 

*i = J 0 < 10 > 

Results of the integrations were used to determine the average fuel-air ratio and the size 
of the undisturbed streamtube upstream of the injector that contains the same air mass 
flow rate as the mixing region. The undisturbed streamtube area Aj was calculated 
for each x/d and q r by assuming that it was of the same shape as A Q ; that is, the 
ratios of the height z to the average width y = A/z of the undisturbed streamtube and 
the mixing region cross section were assumed to be equal. The height of the undisturbed 
streamtube zj was calculated from the continuity equation and the definition of boundary- 
layer displacement thickness: 



Values of zj obtained from equation (11) for the integrated air mass flows from equa- 
tion (10) are given in the tables of figures 12 and 13. In application to the design of 
supersonic combustors, the size and shape of the undisturbed streamtube provide 
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information as to the spacing of injectors and the injector size to obtain penetration to 
the combustor centerline. If it is assumed that injection is from opposite walls and that 
a stoichiometric value of f is desired, the combustor entrance must have a half-height 
equal to the value of zj that corresponds to the x/d station at which a stoichiometric 
average fuel-air ratio was obtained. Confining the flow field by the addition of an opposite 
wall or by the addition of adjacent injectors would be expected to change the mixing rate 
and penetration from that presented herein but would yield a stoichiometric average fuel- 
air ratio. The injector spacing is the average width of the undisturbed streamtube y^. 

An estimation of these parameters is discussed in the next section. 

Estimation of Cold Flow Mixing Parameters 

The mixing length relative to the combustor entrance height and the injector spacing 
required to give a stoichiometric average fuel-air ratio (f = 0.0293 for hydrogen in air) 
have been estimated by superimposing the flow-field cross sections of single injectors. 

A schematic for a two-dimensional configuration with staggered injection from both walls 
is presented in figure 20. Injector spacing was selected so that the half-maximum con- 
centration contour of the single jet flow field was approximately coincident with the half- 
maximum concentration contour of the adjacent and opposite injectors. With this choice 
of spacing the air mass flow that mixes with the hydrogen from one injector is contained 
within the half-a max contour, since the area bounded by the adjacent contours covers the 
entire cross section of the combustor. By denoting the area contained by the half-a max 
contour as A5, equation (10) for the air mass flow becomes 

*l-<W A 5j' 0 <12) 

The height of the confined air streamtube is given by equation (11) by replacing Aq and 
z 0 with A5 and Z5, respectively. 

Values of average fuel-air ratio obtained from the total injected hydrogen mass 
flow and the integrated air mass flow for the unconfined jet (eq. (11)), and the ^/^rnax 
of 0.5 contour (eq. (12)) are presented in figure 21 as a function of x/dj and q r . The 
effect of q r on the decay of f for both conditions was found to be such that a factor 
of q"*/^ provided a reasonable correlation. For the unconfined jet, the value of f 0 
is below stoichiometric for all x/d stations considered. Downstream of 30 injector 
diameters, the decay of f 0 is inversely proportional to approximately x/dj * . For 
the case of simulated opposite wall injection, the average fuel-air ratio £5 is consider- 
ably higher than f 0 and becomes stoichiometric at an x/dj of about 200 for a q r of 
1.0. It is apparent that the lowest value of q r will provide a stoichiometric average 
fuel-air ratio in the shortest distance. Also, using a conincident concentration contour 
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with a value less than half maximum would yield a stoichiometric value of f at a 
shorter x/d since the resulting curve for the decay of f would lie somewhere 
between the curves for f Q and F 5 . However, the uniformity of the flow field would 
need to be investigated to select the optimum value. 

To determine the combustor size, the average fuel-air ratio of the simulated com- 
bustor flow field is plotted in figure 22 as a function of x/zj. The value of zj repre- 
sents the combustor entrance half-height. From figure 22, £5 is stoichiometric at a 
value of x/zj of approximately 40. The hydraulic diameter of a two-dimensional duct 
is twice the duct height or 4zj. The mixing length required for a stoichiometric aver- 
age fuel-air ratio is approximately 10 hydraulic diameters which agrees with rule-of- 
thumb values sometimes considered for engine design. For a q r of 0.5, the value of 
zj for a stoichiometric value of £5 is approximately 3.3 injector diameters — corre- 
sponding to an x/dj of 140 from figure 21. The value of yj , which is the required 
injector spacing, corresponding to these conditions is approximately 3.0 injector diam- 
eters. The analysis just described is of interest in illustrating trends and general mag- 
nitudes; however, the numerical data values cannot be expected to apply directly to 
scramjet combustor design. 


CONCLUDING REMARKS 

An investigation of the penetration and mixing of hydrogen gas injected normal to a 
Mach 4.03 air stream from a sonic injector has been conducted for ratios of jet dynamic 
pressure to free-stream dynamic pressure q r of 0.5 to 1.5. A correlation provided 
the best agreement with the present data (which correlated as a function of the 0.3 power 
of q r ^ whereas other correlations underpredicted the effect of q r . Differences between 
the present data and data correlations in the literature are probably due, in part, to the 
relatively thick boundary layer used in this investigation. 

The decay of the maximum concentration, correlated with the potential core length 
parameter, was inversely proportional to approximately the square root of downstream 
distance at locations less than 50 injector diameters. Values of the potential core length 
agreed well with a correlation which gave the potential core proportional to the 0. 54 power 
of the ratio of jet mass flux to free-stream mass flux. 

Examination of the nondimensional concentration profiles on the vertical center line 
suggest that the profile shape above the point of maximum concentration is not dependent 
on the dynamic -pressure ratio and may be represented by a Gaussian -type function. The 
profiles show similarity at values of the downstream distance x/d equal to or greater 
than 60. Horizontal concentration profiles through the point of maximum concentration 
are also represented by a Gaussian-type function and exhibit similarity at values of x/d 
less than 60. 
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The loss in momentum and total pressure of the airstream caused by turning and 
accelerating the hydrogen jet resulted in an extensive region of very low total pressure. 
The region extended over 60 percent of the mixing region height near the center of the 
flow field and 30 percent of the mixing region width with a total pressure less than 8 per- 
cent of the free stream. 

The results obtained from the present data were used to simulate staggered 
opposite -wall injection in a two-dimensional supersonic combustor by superimposing 
the single jet flow fields. An estimation of the cold flow mixing parameters indicated 
that the mixing length required for a stoichiometric average fuel -air ratio was approx- 
imately 10 hydraulic diameters. 

Langley Research Center. 

National Aeronautics . ; ce Administration, 

Hampton, Va., January 5, 1971. 
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Figure 1.- Sketch of model. All linear dimensions are in centimeters 
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Figure 2.- Boundary -layer surveys on flat -plate model at injector exit station. No injection; T = 300 K. 









Figure 3.- Tunnel and injector flow schematic. 





linear dimensions are in millimeters. 
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(a) Maximum concentration trajectory. 



(x/d) q r " -8 


(b) Half -maximum concentration trajectory. 



(c) Penetration trajectory; v - 0.005. 

Figure 8.- Effect of dynamic-pressure ratio on penetration trajectories 
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(a) Effect of dynamic -pressure ratio. 

Figure 9.- Decay of maximum concentration with downstream distance. 
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Figure 9.- Concluded. 
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Figure 10.- Variation of potential core -length correlation parameter with mass flux ratio. 
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Figure 11.- Nondimensional concentration profiles. Vertical survey. 
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(c) q r = 1.5. 
Figure 12.- Concluded. 
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Figure 13.- Nondimensional total -pressure profiles. Vertical 
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(a) x/d = 7, 30, and 60. 


Figure 14.- Nondimensional concentration profiles. Horizontal 
survey through point of maximum concentration. 
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(b) x/d = 120 and 200. 
Figure 14.- Concluded. 
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y/d - (y/d) a 
(y“d ) 0 


(c) q r = 1.5. 

Figure 15.- Nondimensional velocity profiles. Horizontal survey through point 

of maximum concentration. 






y/d - <y/d) a 

(yTd ) 0 

(b) q r = 1.0. 

Figure 16.- Nondimensional total -pressure profiles. Horizontal 
survey through point of maximum concentration. 
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18.- Continued. 
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Figure 21.- Average fuel-air ratio decay. 
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Figure 22.- Decay of average fuel -air ratio for simulated opposite wall injection. 
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